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Abstract 

We investigate ultraviolet fixed points in the real-time evolution of non-Abelian gauge fields. 
Classical-statistical lattice simulations reveal equal-time correlation functions with a spectral 
index k = 1.5. Analytical understanding of this result is achieved by employing a 2PI- loop 
expansion for the quantum theory. 
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1. Introduction 

The occurrence of turbulent power law spectra after the saturation of chromo-Weibel 
instabilities has been established by several authors [1]. However, the value of the spec- 
tral index characterising these spectra is still under debate. 

We investigate this question in two ways [2]: Firstly, we carry out a perturbative, an- 
alytical calculation for the UV- part of the spectrum similarly as in [3l. This yields a 
spectral index k = 1.5 in the equal-time correlation function of the gauge fields. Subse- 
quently, we run classical-statistical numerical simulations starting from nonequilibrium 
initial conditions. We find that the spectral index from the numerical solution agrees 
with the result from the analytical computation. 
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2. Analytics 



Our interest in this work focuses on the anti-commutator expectation value 

F^ti^,y)^^{{A;{x),Atiy)}) (1) 

of the gauge field. It appears in a decomposition of the full propagator according to 

G;1{x, y) ^ F;t{x, y) - '-pf^ix, y) sign(x° - y") (2) 

where plfl{x,y) = i( [A" (a;), ^^(y)] ) is the spectral function. The nonlocal part of the 
self-energy can be decomposed in a similar way as the propagator [4], 

niiix, y) ^ ^,{x, y) - ^n("^5^,(a;, y) sign(x° - 2/°) , (3) 

with (anti-)symmetric H^^^ ( 11^^^ ). All these objects are diagonal in the colour indices 
and thus we will drop the latter in what follows. 

The goal of this work is to study fixed-point solutions. These are invariant under 
time translations and we additionally assume spatial translation invariance as well. This 
implies that all correlation functions and self-energies only depend on the displacement 
X — y, that is F^^{x,y) = Ff^,j{x — y) = Fp^v{x^ — y^,x — y), and analogously for the 
other objects. A solution with these properties satisfies the following identity in Fourier 
space [5]: 

Tl^^]{p)F,.{p) - 11^7) iP)PyAp) = . (4) 

In order to find an analytical solution we will later integrate this equation over p. 

Now we make the following scaling ansatz for the symmetric and anti-symmetric cor- 
relation functions: 

Fp^{sp) = and Pt^t.{sp) = s~2pp^(p) . (5) 

The partial Fourier transform of such a solution for F possesses the property 

F^,(A< = 0,sp) = s-(i+'')F^,(At = 0,p) . (6) 

In a perturbative calculation the dominant contribution to the self-energies in the UV 
stems from the diagram depicted in Fig. [T] Using the perturbative three- vertex for the 
SU(2)- gauge theory we can write the self-energies as 

n(7)b) = V/ / ^Sip + k + q){F^pik)F,4q)-p^pik)p,4q)} 

X [g^^ip - k)^ + g^Hk - qf + g^^^{q - pf ] [g'^'ip - k f + g^^k - qY + g"^ [q - pf ] 

n(^)(p) = V/ / -^J{P + k + q){F.p{k)pU<l)+P.p{k)F^A<l)] 

X [g^^{p - kf + g-\k - qf + g^^q - pf ] [g^'iP ^ W + - qY + g"'' [q - pf] . 

(7) 

In the classical limit where F'^ ^ p^ the contribution to ^l^p-^ of the two spectral functions 
p in the curly brackets vanishes. Inserting the above expressions for the self-energies 

2 



Fig. 1. The leading order contribution to the gluon self-energy in a perturbative expansion. The lines 
represent full propagators. 

into Eq. ([4]), integrating over the spatial p and carrying out a generalised Zakharov 
transformation similarly as in [3] we obtain 



Further details of the calculation can be found in [2] . We emphasise that without the 
classicality assumption there would not be a fixed-point solution. 

3. Numerical results 

Our numerical setup is the same as in [6]. We briefly summarise the details that are 
relevant for the present discussion. 

We use a standard lattice discretisation scheme in terms of link variables Ux,i_i repre- 
senting parallel transporters from the lattice site x to x jl. Given the link variables one 
can compute the gauge fields as x) ~ —i Tr [(j''U^{t, x)) / {2ag) where cr'' , 6 = 1, 2, 3, 
are the Pauli matrices and a is the lattice spacing. 

We compute the real time evolution of gauge fields by solving the discrete equations of 
motion derived from a variation of the Wilson lattice action in Minkowski spacetime. The 
initial conditions are sampled from a Gaussian distribution such that { \ A^{to,p) ) oc 
exp{— 0.5(p^ +Py)/A'^ — 0.5(p^)/A^} for the gauge fields in momentum space. The ex- 
pectation values occurring in the definitions of the correlation functions ([1]) have to be 
taken with respect to this distribution. 

We numerically follow the time evolution of the gauge fields to times far beyond the 
saturation of the instabilities. Computing the Fourier spectrun{3 of Ffj_^{At — 0,x) — 
{ Afj_{t, x)A^{t, 0) } at late times in Coulomb gaugj^ we observe in Fig. [2] that the spec- 
trum converges towards a stationary curve. Fitting a power law ansatz F^i/(0,p) = 
C|p|~*-^+"-' as in Eq. ((6]) to the latest available spectrum on the domain 0.6 < p- e^^l^ < 
2.0 we obtain 

K = 1.52 ±0.07. (9) 

We can produce a slightly higher exponent of k = 1.69 ±0.03 by including also the lowest 
momenta into the fit. However, this increases per degree of freedom to 1.7, in contrast 
to x^/ d. o. f. =1.1 for the fit underlying Therefore, the result quoted in ^ is more 
reliable and we note that it is in excellent agreement with Eq. ([5]). In any case our data 
is incompatible with k = 2 and also with the thermal value k = \. The fit results are 
visualised in Fig. O 



^ Precisely speaking, _F^^ here denotes the classical approximation to the object defined in Eq. jlj. 
^ Our gauge fixing employs the stochastic overrelaxation algorithm described in [?]• 
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Fig. 2. Spectra of Ffiv{A.t = 0, to) = {Afi{t,x)Av{t,0)) , fj, = u . Grey symbols: Spectra at times 
t = 212 - 705 £-1/4. Black squares: t = 940 e ^Z'*. The blue lines represent the fit results, confer the text 
for details. The momenta are perpendicular to the spatial direction of the gauge fields. 

4. Conclusions 



The analytical and numerical results for the spectral index are in excellent agreement. 
Their numerical value also agrees with the one found for the UV- spectrum of scalar field 
theories in the context of early-universe cosmology [318] . 

What is the meaning of these findings? Comparing the present results to IS'S' shows 
that some features of QCD are quantitatively the same as for scalar theories. This might 
allow to establish the notion of universality far from equilibrium. In fact, quantum cor- 
rections will eventually drive the system away from the classical fixed-point solution in 
the UV. However, they will not alter the infrared behaviour, as it is known from scalar 
theories [3]. It remains as a future project to search for nonthermal fixed points in the 
IR- sector of the gauge theory. Also, it would be desirable to discuss fixed points in terms 
of gauge-invariant quantities. 

I thank my collaborators Jiirgen Berges and Denes Sexty for their contributions to this 
project. 
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